function input_image = arnold_inverse_transform(output_image, iterations)
    % Arnold逆变换
    [h, w]=size(output_image) ;
    input_image = output_image;

    a = 3;
    b = 4;
    for i=1:iterations
        for y=1:h
            for x=1:w            
                xx=mod((a*b+1)*(x-1)-b*(y-1),h)+1;
                yy=mod(-a*(x-1)+(y-1),h)+1  ;        
                input_image(yy,xx)=output_image(y,x);                   
            end
        end
        output_image=input_image;
    end

end